Threshold-Resummed Cross Section for the Drell-Yan Process 



in Pion-Nucleon Collisions at COMPASS 
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We present a study of the Drell-Yan process in pion-proton collisions including next-to-leading- 
logarithmic threshold-resummed contributions. We analyze rapidity-integrated as well as rapidity- 
differential cross sections in the kinematic regime relevant for the COMPASS fixed target experiment. 
We find that resummation leads to a significant enhancement of the cross section compared to fixed- 
order calculations in this regime. Particularly large corrections arise at large forward and backward 
rapidities of the lepton pair. We also study the scale dependence of the cross section and find it to 
be substantially reduced by threshold resummation. 

PACS numbers: 12.38.Cy,12.38.Bx,14.40.Be 



I. INTRODUCTION 

Drell-Yan lepton pair production by incident charged 
pions is still the only available hadronic process that al- 
lows for the extraction of the internal "partonic" struc- 
ture of the pion in the valence region. Unfortunately, ex- 
perimental data for this process are rather scarce [l], ■ 
Several analyses of the available data have been per- 
formed using leading order (LO) or next-to-leading order 
(NLO) partonic cross sections, with the goal of extract- 
ing the parton distribution functions of the pion ■ 
The resulting LO or NLO valence distributions showed a 
linear (~ (1 — x) 1 ) or slightly faster fall-off at high x, at 
odds with theoretical predictions based on perturbativc 
QCD m , and calculations using Dyson- Schwinger equa- 
tions [3], which prefer a much softer behavior ~ (1 — x) 2 
at high x. For a review, see Q. 

In a recent analysis @ we showed that a valence distri- 
bution with a fall-off ~ (1 — x) 2±01 is in fact well consis- 
tent with the Drell-Yan data, if large logarithmic contri- 
butions arising near the threshold for the partonic reac- 
tion are taken into account to all orders in perturbation 
theory. These so-called "threshold logarithms" strongly 
enhance the partonic cross section near threshold, so that 
a softer valence distribution is sufficient to describe the 
data. We also found in Ref. [9( that the available pionic 
Drell-Yan data are not able to completely determine the 
valence parton distribution of the pion. The data are 
almost equally well described by a valence distribution 
that carries 60 %, 65 % or 70 % of the pion's momentum 
at a low input scale Qo = 0.63 GeV, as long as it behaves 
approximately as (1 — x) 2 at high x. The quality of three 
corresponding fits that we performed in |9| only differed 
by about one unit in x 2 - This uncertainty in the valence 
momentum is also manifest in the earlier NLO analy- 
ses. At Q = 2 GeV the valence parton distribution of 
SMRS carries 46 % of the pion's momentum, whereas 
the valence distribution of GRS Q| carries only 40 %, 
although both distributions describe the same data sets 



equally well. The origins of this ambiguity are the large 
overall systematic uncertainties of the Drell-Yan data. 
These uncertainties can be best discussed by introduc- 
ing an additional if-factor that multiplies the theoretical 
cross section. The numerical value of this if-factor is 
strongly correlated with the first moment of the chosen 
valence distribution (see Table I of Ref. [9J). Hence, pion 
Drell-Yan data with a well understood normalization are 
urgently needed to really pin down the pion's valence 
distribution. 

The up coming fixed-target ttN Drell-Yan experi- 
ment [Of at COMPASS is hoped to resolve this is- 
sue. In this paper, we present detailed predictions for 
the rapidity-integrated as well as the rapidity-differential 
cross section for the kinematics relevant at COMPASS. In 
the light of our study , it is expected that threshold log- 
arithms will also play a significant role for the Drell-Yan 
cross section at COMPASS and lead to large corrections. 
We will therefore base our predictions on the thresh- 
old resummation technique. We note that there have 
been numerous earlier phenomenological applications of 
threshold resummation in the Drell-Yan process 
both for fixed-target and for collider energies. The spe- 
cific application to ttN scattering, and the resulting phe- 
nomenological studies for COMPASS, are the new ele- 
ments of this paper. 

The remainder of this paper is organized as follows. 
In Sec. HI] the basic framework for the calculation of 
the rapidity-differential Drell-Yan cross section is pre- 
sented. We discuss fixed-order corrections as well as the 
full next-to-leading logarithmic resummation of thresh- 
old logarithms. This includes the discussion of Mellin 
and Fourier moments of the cross section, which are 
useful tools for threshold resummation. In Sec. IIIII 
we present phenomenological results for the kinematic 
regime of the Drell-Yan experiment at COMPASS. We 
show both rapidity-integrated and rapidity-differential 
cross sections. We finally draw our conclusions in Sec. lIVI 



II. THEORETICAL FRAMEWORK 

We follow in this section the framework laid out in 
Ref. [Hj], where threshold resummation effects in W^- 
boson production at hadron colliders were studied. We 
consider the inclusive cross section for the production of 
a /x + /Lt _ pair of invariant mass Q and rapidity 77 in the 
process 

n-(P 1 )p(P 2 )^ fx+fx-X, (1) 

where Pi and P2 are the four-momenta of the initial- 
state particles. According to the relevant factorization 
theorem, at high Q the rapidity-differential cross section 
may be written as 

da f 1 d Xl f 1 dx 2 „, 2 \tp, 2n 

dQ 2 dr] ~^Jx° x i Jx% x z 

Xe a b^ab(xi,Xi,X2,X2,Q/fJ-), (2) 

where ct = 4ira 2 /9Q 2 S with S = {P1+P2) 2 the hadronic 
center-of-mass energy squared. The coupling e a b equals 
e 2 for the qq and qg, qg scattering processes which we 
are interested in, where e q denotes the quark's fractional 
electromagnetic charge. In terms of the rapidity rj the 
lower bounds of the x\ and X2 integrals are 

<a = v^e^, (3) 

with r = Q 2 /S. The sum in Eq. © runs over all par- 
tonic channels, with /J and f? the corresponding parton 
distribution functions (PDFs) of the pion and the pro- 
ton, and uj a i, the hard-scattering functions. The latter 
can be computed in perturbation theory as series in the 
strong coupling constant a s . The parton distribution 
functions as well as the hard-scattering functions depend 
on the factorization and renormalization scales, which we 
choose to be equal and collectively denote as \i. 

At leading order O{a s ) only the quark- antiquark anni- 
hilation channel qq — > 7* — > contributes, for which 
one has in our normalization 

uf-j = xix 2 S(xi - x\)8(x 2 - x 2 ). (4) 

At next-to-leading order, apart from the 0(a s ) correc- 
tions to the qq process, also additional processes con- 
tribute to the cross section, namely qg — » j*q and 
qg — > 7*9. The NLO partonic cross sections in the MS 
scheme, which is the scheme we adopt throughout this 
work, can be obtained from [3. They are also collected 
in the Appendix of Ref. [3]. 1 



We note that Refs. [1 fill adopt a non-standard polarization av- 
erage for incoming gluons in dimensional regularization. In order 
to correct for this, one simply needs to multiply the arguments 
of the logarithms in Eqs. (A8) and (A20) of Q by e _1 fist . 



As mentioned in the Introduction, the Drell-Yan cross 
section receives large logarithmic corrections near the 
threshold for the partonic reaction [2(| . This threshold is 
defined by z — Q 2 /x\x 2 S = 1, where x\ and X2 are the 
momentum fractions of the partons participating in the 
hard-scattering reaction. As z increases towards unity, 
most of the initial partonic energy is used to produce 
the virtual photon. Therefore, little phase space remains 
for real-gluon radiation, while virtual-gluon diagrams 
may still contribute fully. The infrared cancellations be- 
tween the virtual and the "inhibited" real-emission di- 
agrams then leave behind large logarithmic corrections 
to oj q q. The leading terms among the resulting thresh- 
old logarithms are of the form a^ln 2fc-1 (l — z)/(l — z) 
at the fcth order of perturbation theory. Subleading 
terms are down by one or more powers of the logarithm. 
The threshold logarithms become particularly important 
when r = Q 2 /S, the hadronic analog of z, is large, which 
is generally the case in the fixed-target regime. The fact 
that the parton distribution functions are steeply falling 
functions of X\ or X2 emphasizes the threshold region 
in the cross section even for values of t substantially 
smaller than one. In this kinematic regime the logarithms 
ln 2fc ~ 1 (l — z)/(l — z) compensate the smallness of a k sl and 
it becomes necessary to resum the large corrections to all 
orders in the strong coupling. Such "threshold resumma- 
tion" has originally been derived for the Drell-Yan pro- 
cess and deep inelastic scattering a long time ago |20| . 
The techniques developed in these seminal papers have 
been extended and successfully applied to the resumma- 
tion of large logarithmic contributions in numerous other 
hard QCD processes. 

Threshold resummation may be achieved in Mellin mo- 
ment space, where phase space integrals for multiple- 
soft-gluon emission decouple. For the rapidity depen- 
dent cross section, it is convenient to also apply a Fourier 
transform in 77 [13, [21| (alternatively, one can also use a 
double Mellin transform [22[ ). Under combined Fourier 
and Mellin transforms of the cross section, 

„{N t M)= fdrr^ /' n * d^-g-, (5) 
Jo J- in -i. dQ 2 dr] 

the convolution integrals in © decouple into ordinary 
products [IHHlj]. Defining the moments of the PDFs, 

f N (n 2 )= [ dxx N - 1 f{x,^), (6) 
Jo 

and introducing the corresponding double transform of 
the partonic hard-scattering cross sections, 

Q ab {N,M)= [ 1 dzz N - 1 ^ drK Mj W, (7) 
Jo J- In 4= 

where fj = rj — \ 111(0:1/0:2) is the partonic center-of-mass 
rapidity, one finds: 

a{N,M) = a J2fa N+i " e ab u ab {N,M). (8) 

a , b 



The double transform thus factorizes the PDFs and the 
perturbatively calculable hard-scattering functions. The 
leading order contribution to the hard-scattering func- 
tion £j a b(N, M) is easily calculated by making use of the 
relations 



^ = sfle 

%2 



We obtain from Eq. ^ for the Fourier transform of i 



(9) 



,(o). 



ln(l/Vi) 

dfje^x^Sixi - x° 1 )S(x 2 - x° 2 ) 

-hx(l/VF) 

In(VV^) (10) 

drje iMr >S(l - ^e v )S(l - Vie - ") 

•MW5) 
= cos (M ln(l/v/i)) 5(1 - z) . 

Here we have appropriately averaged over the two possi- 
ble solutions for the integral. The emerging factor 8(1— z) 
in Eq. (JTUJ) is just the LO hard-scattering contribution 
to the rapidity- integrated Drell-Yan cross section. Hence, 
the Fourier transform of the LO rapidity-differential par- 
tonic cross section is equal to the LO rapidity-integrated 
partonic cross section times cos (Mln(l/y / i)) [12|]. In 
the near-threshold limit z — > 1 this cosine factor becomes 
subleading: 

cos (Af ln(l/Vi)) = 1 - f + 0((1 - z) 4 M 4 ). 

8 

(11) 

As was discussed in Refs. [T3, E3, HH, even at 
higher orders the dependence of the double moments 
u a b=qq{N, M) on M becomes subleading near thresh- 
old, whereas the iV-dependence is identical to that of 
the rapidity-integrated cross section. Therefore, the re- 
summed expression for u a b =qq (N, M) is equal to that 
for the total (rapidity-integrated) cross section. It was 
shown in Rcf. [12J that keeping the cosine term in 
Eq. (fl0|) in the resummed Co a b= qq slightly more faithfully 
reproduces the rapidity dependence at each order of per- 
turbation theory. This can be easily achieved by writing 
the cosine as 



cos(Mln(l/VS)) = \ 



JM/2 



-iM/2 



(12) 



which, when combined with Eq. ([7]), leads to a sum of 
two terms with Mellin moments shifted to N ± iM/2. 

Threshold resummation for the Drell-Yan process re- 
sults in the exponentiation of the soft-gluon correc- 
tions. To next-to-leading-logarithmic (NLL) order the 
resummed cross section is given in the MS scheme by 



i-C 



Inw, 9 = C q ( ^,a s ( M 2 )J + 2 J dC 

Ki-C) 2 Q 2 dk 2 
" 1 -j^A q (a s (kx)), 



(13) 



where A q (a s ) is a perturbative function, the O(ai ) part 
of which is sufficient for resummation to NLL 1201: 



(14) 



with 124 



q 



q 



:C f 



C A 



67 
18 



6 



9 



-N 



(15) 

Here C F = 4/3, C A = 3. The first term in Eq. (13) does 
not originate from soft-gluon emission but instead mostly 
contains hard virtual corrections. It is also a perturbative 
series in a s , and we only need its first-order term: 



C a 



-C P -4 



2f_ 
3 



3, Q 
-In — 

2 /i 



+ 0(a 2 s ), (16) 



whose exponentiated form has been established in Ref. 

Since the perturbative running coupling a s (k±) di- 
verges at k± — Aqcd, Eq. (fT3) as it stands is ill-defined. 
The perturbative expansion of the expression shows fac- 
torial divergence, which in QCD corresponds to a pow- 
erlike ambiguity of the series [26]. It turns out, however, 
that the factorial divergence appears only at nonleading 
powers of the momentum transfer. The large logarithms 
we are resumming arise in the region (2Cj z < 1 — 1 / N in 
the integrand of the second term in Eq. (fT5|) . One there- 
fore finds that to NLL they are contained in the simpler 
expression 



- ; dk\ 

'■ /N 2 "<_L 

dk 2 



-g±A q (a s (k ± ))ln 



Nk, 



j2 A q (a s (k ± ))lnN 



(17) 



for the second term in Eq. (fT3|) . where N — Ne lB with 
the Euler constant je- This form is used for "mini- 
mal" expansions [27| of the resummed exponent. From 
Eq. p7|) one obtains for the resummed exponent to NLL 
accuracy [27l[28j: 



In ui qq = C q + 2/i (1) (A) In N + 2h {2} (x, 



where 



A = b a s {^ 2 )\nN. 



(18) 



(19) 



The functions collect all leading- logarithmic 

and NLL terms in the exponent, which are of the form 



a k s ln fe+1 N and a k s ln K N, respectively. They read 



ft( 2 ) A 



Q 2 



.4, 



(i) 



2tt6 A 



[2A + (1 - 2A) ln(l - 2A)] 



(2) 



^[2A + ln(l-2A)] 



27T&3 



[2A + ln(l - 2A) 



ln 2 (l - 2 A) 



(i) 



+^r[^ + Ml-2A)]ln^ 



2iTb 



\n(N)\n%, (20) 



The last term of the function depends on the fac- 
torization scale and compensates the evolution of the 
parton distribution functions. The scale dependence of 
the second-to-last term results from the running of the 
strong coupling constant. Since scale evolution exponen- 
tiates and is therefore taken into account to all orders, 
one expects [1^, [3(| a significant decrease in the scale de- 
pendence of the resummed cross section compared to a 
fixed order cross section. 



where 



b 



— (11C A - 2N f ) (21) 

1Z7T 

(l7C A - 5C A N f - 3C F N f ) . (22) 



24tt ; 



As was shown in Refs. [3ll - l33| . it is possible to im- 
prove the above formula by taking into account certain 
subleading terms in the resummation. As in Ref. [l2T | we 
rewrite Eqs. CED-flU]) as 



hlLO a 



1 / 4 (1) 

- [2A + ln(l-2A)] — ^ 



A q % 




-C F 



2-k 2 



27r6g 



(23) 



{-2A^\nN-B^} (^l\n% 



ln(l - 2 A) 

7T&Q 



r 



where Bq = 



3C F /2. The last term in Eq. ([23]) is the 
leading-logarithmic expansion of the integral 



L fc 2 7T [ 



N - 



(24) 



The term in square brackets is the leading term in the 
large- N limit of the anomalous dimension of the one-loop 
diagonal (q q) splitting function P^, i.e. it governs 
the evolution of the parton distributions between scales 
H and Q/N. Replacing it by the full flavor nonsinglet LO 
splitting function (32j |. 



In AT 



C F 



25i (A) 



1 



A(A + 1) 



(25) 

fully reproduces the diagonal part of the quark and anti- 
quark evolution. This replacement further reduces the 
scale dependence of the resummed cross section. We 
could also include a non-diagonal contribution from g — > 
q splitting, corresponding to singlet mixing. However, 



this contribution turns out to be numerically unimpor- 
tant for the Drell-Yan process in the present kinematics. 

As the exponentiation of soft-gluon corrections is 
achieved in Mellin moment and Fourier space, the 
hadronic cross section differential in Q 2 and rj is obtained 
by taking the inverse Mellin and Fourier transforms of 
Eq. ©: 



da 



dQ 2 d7y 



AM 
~2w ( 



-iMrj 



C+ioo 



C—ioc 



2~ki 



-N 



o(N,M). 



(26) 

When performing the inverse Mellin transform, the pa- 
rameter C usually has to be chosen in such a way that all 
singularities of the integrand lie to the left of the integra- 
tion contour. The resummed cross section, however, has 
a Landau singularity at A = 1/2 or N = exp(l/2a s &o), 
as a result of the divergence of the running coupling a s 
in Eq. ([13]) for k± — > Aqcd- For the Mellin inversion, we 
adopt the minimal prescription developed in Ref. [13] 
to deal with the Landau pole. For this prescription the 
contour is chosen to lie to the left of the Landau sin- 



gularity. Above and below the real axis, the contour is 
tilted into the half-plane with negative real part. This 
improves the convergence of the integration, since con- 
tributions with negative real part are exponentially sup- 
pressed by the factor t~ n in Eq. (|26|) . As mentioned 
earlier, the moment-space singularities of the parton dis- 
tribution functions are shifted by ±iM/2 from the real 
axis due to the Fourier transform. Ref. (2l| provides de- 
tails for how to prevent the tilted contour from passing 
through or below those singularities. We note that an 
alternative possibility for dealing with the Landau singu- 
larity is to perform the resummation directly in z-space 

We match the resummed cross section to the NLO one 
by subtracting the 0(a s ) expansion of the resummed ex- 
pression and adding the full NLO cross section [l2[. This 
"matched" cross section consequently not only resums 
the large threshold logarithms to all orders, but also con- 
tains the full NLO results for the qq and qg channels. We 
will occasionally also consider a resummed cross section 
that has not been matched to the NLO one. We will refer 
to such a cross section as "unmatched" . 



III. PHENOMENOLOGICAL RESULTS 

We now present our numerical results for the Drell-Yan 
cross section at COMPASS. We will consider both the 
rapidity-integrated and the rapidity-differential hadronic 
cross section. Our main goal is to investigate the size of 
the threshold resummation effects. The ir~ beam fore- 
seen at COMPASS has an energy of 190 GeV. It is scat- 
tered off a proton target at rest, so that the resulting 
center-of-mass energy of the system is \fS ~ 19 GeV. 

For the pionic parton distribution functions we use the 
ones for the "preferred fit" of our previous study Q- 
We remind the reader that these were extracted from a 
fit to the earlier 7riV Drell-Yan data 0, 0], using NLL 
threshold-resummed cross sections in the MS scheme. 
For the proton target we use the NLO (MS scheme) 
CTEQ6M [H parton distributions. Unless stated oth- 
erwise, we choose the renormalization and factorization 
scales as [i = Q. 

We start by considering the differential cross section 
da/dQ, integrated over all rapidities, to show the rele- 
vance and the validity of the resummation effects over 
the whole range of the invariant mass Q. Here we ignore 
for simplicity charmonium and bottonium resonances in 
the lepton pair spectrum, whose contributions are dom- 
inant for resonant invariant masses, and calculate only 
the smooth (continuum) part of the cross section. Fig- 
ure Q] shows the cross section Q 3 da/dQ at \/~S = 19 GeV 
at fixed order (LO and NLO), as well as for the NLL- 
resummed case. It can be seen that resummation leads 
to a significant enhancement of the cross section over LO, 
which increases strongly with invariant mass. This be- 
comes even more apparent in Fig. [21 where we show the 
"K- factor" , defined as the ratio of the cross section to the 




10 

Q (GeV) 

FIG. 1: Rapidity-integrated Drell-Yan cross section Q 3 da/dQ 
for 7T~ p scattering at \/S = 19 GeV, at LO, NLO and NLL- 
resummed, as a function of the invariant mass Q of the lepton 
pair. 



LO one: 



K = 



da/dQ 
d<7 LO /dQ 



(27) 



The "K-factor" is plotted for the NLO and the NLL- 
resummed result. We also expand the unmatched re- 
summed cross section in powers of a s . The results for the 
first, second and third order expansion are also shown in 
Fig. [2J One can see that in the fixed-target regime higher 
orders (beyond NLO) still make large contributions to the 
cross section, especially at high invariant mass Q. This 
finding is in line with that in the earlier study 11[ for pp- 
scattering. We also observe that the exact NLO cross sec- 
tion agrees extremely well with the first order expansion 
of the unmatched resummed result. This demonstrates 
that the logarithmic contributions from soft gluon radi- 
ation, which we resum to all orders, give by far the most 
important contribution to the cross section, not only very 
close to threshold as r = Q 2 /S —> 1, but also for rather 
moderate values of r. 

Next, we present the results for the rapidity distri- 
butions da /dQdrj. As mentioned above, charmonium 
and bottonium resonances complicate the calculation of 
Drell-Yan cross sections. Therefore usually only lepton 
pairs with invariant mass Q between the J/\& and T 
resonances and above the T are considered. Since the 
Drell-Yan event rate decreases rapidly with y/r, it may 
not be possible to measure it above the T resonance in 
the medium-energy fixed-target regime accessed by the 
COMPASS experiment. We therefore make predictions 
for y/r — 0.3 and yfr = 0.45, corresponding to Q = 5.7 
GeV and Q — 8.6 GeV, respectively. Our results are 
presented in Figs. [3] and [H Again the resummed cross 
section and the fixed-order NLO and LO ones are shown. 
As before, we expand the unmatched resummed result 
in powers of a s and find that the first order expansion 
agrees very well with the exact NLO result for yfr — 0.45. 




FIG. 2: "K-factors" as defined in Eq. IpTJl at </S = 19 GeV 
as functions of the lepton pair mass Q, at NLO (symbols) and 
for the NLL-resummed case. Also shown are the expansions 
of the resummed cross section to first, second and third order 
in the strong coupling. 



50 



9 40 

O 

a 30 

8*20 



NLL resummed 
1st order expansion 
NLO 
LO 



10 




FIG. 3: Rapidity-differential Drell-Yan cross section 
dcr/dA/d?7 for ir~ p scattering at VS = 19 GeV and = 
0.3. The LO, NLO and NLL-resummed cross sections as well 
as the first order expansion of the unmatched NLL-resummed 
cross section are shown as functions of the rapidity rj of the 
dimuon pair. 



For yfr = 0.3, further away from threshold, the first or- 
der expansion of the threshold resummed cross section 
lies very slightly below the exact NLO result for central 
rapidities. This is due to the fact that the contributions 
from the threshold region z — > 1 do not entirely dominate 
the cross section in this rapidity regime. As expected 
from our results in Figs. [1] and [2j at fixed rapidity the 
threshold resummation effects become more important 
as r increases, resulting in a fairly large enhancement of 
the resummed cross section at \[r = 0.45. Nevertheless 
significant contributions from threshold resummation are 
still present in the cross section also for relatively modest 
values of r. 



FIG. 4: Same as Fig. O but at yfr = 0.45. 



It is also interesting to examine to what extent resum- 
mation affects the shape of the rapidity dependent cross 
section. Figure [5] shows the ratios 



( do- rcs \ 

/ dcr LO \ ' 



K 



NLO 



/ dffNLO 

V dQd?? 
/ da 



LO 



\dQdrj 



(28) 



as functions of the pair rapidity, at ^/S =19 GeV and 
y/r = 0.3. One can see that K res becomes very large 
towards the boundaries of the r\ interval. The resummed 
cross section shows a particularly sizable enhancement 
above the NLO one at high rapidities. This enhancement 
is due to the fact that at fixed r the limit rj —¥ r/ max cor- 
responds to the limit z — > 1 at parton level. In this limit 
threshold logarithms become large regardless of the value 
of t. Large r and/or high rapidities in the fixed-target 
regime probe high momentum fractions x in the parton 
distribution functions. Including threshold resummation 
in the analysis of parton distribution functions may hence 
have significant effects on their extracted high-cc behav- 
ior. This was examined recently in the context of the 
Drell-Yan process [9j and deep-inelastic lepton scatter- 
ing El. 

The crucial quality test for any higher order calcula- 
tion is the extent to which it reduces the scale ambiguity 
inherent to any perturbative QCD calculation. We ex- 
amine the scale dependences of the rapidity-integrated 
and the rapidity-differential cross sections in Figs. [S] 
and respectively. Again we show the LO, NLO and 
NLL-resummed results at = 19 GeV, now vary- 
ing the renormalization and factorization scales between 
/i = Q/2 and fj, = 2Q. Note that in Fig. [5] we have for 
better visibility multiplied the LO cross section by 1/2 
and the resummed one by 2. Evidently for the integrated 
cross section the scale dependence is decreased by resum- 
mation over the whole range of invariant mass Q, whereas 
going from LO to NLO reduces the scale dependence only 
marginally. Figure [7] shows the scale dependence of the 
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FIG. 5: Ratios K TCS and A"nlo as defined in Eq. ((28)) at yfS = 
19 GeV and yfr = 0.3, as functions of the rapidity r\ of the 
dimuon pair. 




Q (GeV) 

FIG. 6: Scale dependence of the LO, NLO and NLL- 
resummed rapidity-integrated Drell-Yan cross sections at 
\fS = 19 GeV as function of Q. The factorization as well 
as the renormalization scale have been varied between Q/2 
and 2Q. Note that we have multiplied the LO cross section 
by 1/2 and the resummed cross section by 2. 

rapidity distributions at yfr = 0.45. Here we only show 
the NLL-resummed cross section and the NLO one. As 
one can see, the scale dependence is again significantly 
improved by resummation. This applies to all values of 
rapidity; in fact the scale dependence almost vanishes at 
high r] after resummation. 



IV. CONCLUSIONS 
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FIG. 7: Scale dependence of the NLO and NLL-resummed 
rapidity-differential Drell-Yan cross sections at s/S = 19 GeV 
and yfr = 0.45 as function of n. The factorization as well as 
the renormalization scale have been varied between Q/2 and 
2Q. 

resummed cross section to 0(a s ) agrees very well with 
the exact fixed-order calculation. This agreement demon- 
strates that the large threshold logarithms indeed dom- 
inate the Drell-Yan cross section and need to be taken 
into account to all orders. Resumming those logarithms 
leads to a significant enhancement above fixed-order cal- 
culations, even for moderate values of the invariant mass 
Q of the lepton pair. We have also considered the ra- 
pidity dependence of the cross section. We find that 
even in cases where there is only a modest enhancement 
of the rapidity-integrated cross section by resummation, 
the shape of the rapidity- differential cross section is af- 
fected very strongly by resummation at sufficiently large 
forward or backward rapidities. Finally, we have shown 
that the scale dependence of the perturbative cross sec- 
tion is substantially reduced when threshold resummed 
contributions are included. 

Our results overall demonstrate that threshold resum- 
mation effects will be important in the analysis of future 
COMPASS data. While we have only addressed the spin- 
averaged Drell-Yan cross section in this paper, we stress 
that threshold resummation effects are expected to be 
equally relevant also for corresponding spin-dependent 
cross sections, even though they may have a tendency 
to cancel in spin asymmetries. We also note that in the 
light of our study cross sections and spin asymmetries at 
measured transverse momentum q± of the lepton pair, 
which will be a particular focus of the investigations at 
COMPASS, will require additional theoretical consider- 
ation. 



We have presented a phenomcnological study of the 
Drell-Yan cross section for pion scattering off a proton 
target at COMPASS. In the calculation of the cross sec- 
tion we have resummed threshold corrections to next- 
to-leading logarithmic accuracy. The expansion of the 
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